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ABSTRACT 

The ultimate remote sensing benefits of the high resolution Infrared radiance spectrometers will be realized with 
their geostationary satellite implementation in the form of imaging spectrometers. This will enable dynamic 
features of the atmosphere’s thermodynamic fields and pollutant and greenhouse gas constituents to be observed 
for revolutionary improvements in weather forecasts and more accurate air quality and climate predictions. As 
an important step toward realizing this application objective, the Geostationary Imaging Fourier Transform 
Spectrometer (GIFTS) Engineering Demonstration Unit (EDU) was successfully developed under the NASA 
New Millennium Program, 2000-2006. The GIFTS-EDU instrument employs three focal plane arrays (FPAs), 
which gather measurements across the long- wave IR (LWIR), short /mid- wave IR (SMWIR), and visible spectral 
bands. The GIFTS calibration is achieved using internal blackbody calibration references at ambient (260 
K) and hot (286 K) temperatures. In this paper, we introduce a refined calibration technique that utilizes 
Principle Component (PC) analysis to compensate for instrument distortions and artifacts, therefore, enhancing 
the absolute calibration accuracy. This method is applied to data collected during the GIFTS Ground Based 
Measurement (GBM) experiment, together with simultaneous observations by the accurately calibrated AERI 
(Atmospheric Emitted Radiance Interferometer), both simultaneously zenith viewing the sky through the same 
external scene mirror at ten- minute intervals throughout a cloudless day at Logan Utah on September 13, 2006. 
The accurately calibrated GIFTS radiances are produced using the first four PC scores in the GIFTS- AERI 
regression model. Temperature and moisture profiles retrieved from the PC-calibrated GIFTS radiances are 
verified against radiosonde measurements collected throughout the GIFTS sky measurement period. Using the 
GIFTS GBM calibration model, we compute the calibrated radiances from data collected during the moon 
tracking and viewing experiment events. From which, we derive the lunar surface temperature and emissivity 
associated with the moon viewing measurements. 
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1. INTRODUCTION 

The implementation of Infrared ultra-spectral resolution spectrometers on polar orbiting satellites has ushered in 
a new era in satellite remote atmospheric sounding capability. However, the ultimate benefits of the ultra-spectral 
remote sensing technology are expected to be realized through its geostationary satellite implementation in the 
form of imaging spectrometers. This will enable dynamic features of the atmosphere’s thermodynamic fields and 
pollutant and greenhouse gas constituents to be observed for revolutionary improvements in weather forecasts 
and more accurate air quality and climate predictions. As an important step toward realizing this application 
objective, the Geostationary Imaging Fourier Transform Spectrometer (GIFTS) Engineering Demonstration Unit 
(EDU) was successfully developed under the NASA New Millennium Program, 2000-2006. The GIFTS-EDU 
instrument employs three focal plane arrays (FPAs), which gather measurements across the long- wave (685- 
1130 cm -1 ) IR (LWIR), short /mid- wave (1650-2250 cm -1 ) IR (SMWIR), and visible spectral bands. The raw 
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interferogram measurements are radiometrically and spectrally calibrated to produce radiance spectra, which are 
further processed to obtain atmospheric profiles via retrieval algorithms. The absolute radiometric performance of 
the instrument is affected by several factors including the FPA off-axis effect, detector /readout electronics induced 
nonlinearity distortions, and fore-optics offsets. The GIFTS-EDU, being the very first imaging spectrometer 
to use ultra- high speed electronics to readout its large area format focal plane array detectors, operating at 
wavelengths as large as 15 microns, possessed non-linearity’s not easily removable in the initial calibration process. 
In this paper, we introduce a refined calibration technique that utilizes Principle Component (PC) analysis to 
compensate for instrument distortions remaining after the initial calibration procedures. In this paper, we 
emphasize the technique of PC regression calibration and demonstrate the sensitivity of the calibrated radiances 
to atmospheric temperature and moisture variability using data collected with the GIFTS and the AERI, both 
simultaneously zenith viewing the sky through the same external scene mirror at ten-minute intervals throughout 
a cloudless day at Logan Utah on September 13, 2006. In Section 2, we introduce the methodology and results of 
the PC regression calibration technique. In Section 3, we present the temperature and moisture retrieval method 
using the PC calibrated radiances. And in Section 4, we apply the same calibration principle to data collected 
during the moon viewing events. From which, we compute the lunar surface temperature and emissivity. 

2. GIFTS CALIBRATIONS USING PRINCIPLE COMPONENT REGRESSION 

In this section, we introduce a refined calibration technique that utilizes principle component analysis to com- 
pensate for instrument distortions and artifacts remaining after the initial radiometric calibration process. This 
method is applied to data collected with the GIFTS, together with simultaneous observations by the accurately 
calibrated AERI (Atmospheric Emitted Radiance Interferometer), both zenith viewing the sky at ten- minute in- 
tervals throughout a cloudless day of an atmospheric measurement experiment. The PC vectors of the calibrated 
radiance spectra are defined from the AERI observations and regression matrices relating the initial GIFTS 
radiance PC scores to the AERI radiance PC scores are calculated using the least squares inverse method. A 
new set of accurately calibrated GIFTS radiances are produced using the first four PC scores in the regression 
model. The PC-calibrated radiances and initial calibration results are compared to AERI spectra in terms of 
FPA mean and variance radiance differences. 

2.1 GIFTS and AERI Data Definitions 

We define G(n,k,x,y,s) as calibrated GIFTS GBM AVE (Atmospheric Variation Event) long-wave radiance 
spectra, where k = 1, . . . , K is the spectral sample number and K = 778 is the total number of spectral samples. 
The corresponding spectral wavenumber is written as = &k 0 + kAa , where Act = 0.5733 cm -1 is spectral 
resolution and covers from 684.56 cm -1 to 1130.04 cm -1 . The variable n = 1, . . . , N is the temporal sample 
number and N = 52 is the total number of temporal samples; i.e., t = to+nAt, and At = 10 minutes. The variable 
s = 1, . . . , S is the scan number and S' = 5 is the total number of scans collected during each data measurement 
period. Finally, x, y = 1, . . . , 128 represent the pixel index numbers in rows and columns, respectively. We define 
G(n, k , x, y ) as the scan-averaged data cubes, and G(&, x, y ) as the scan-temporal-averaged data cube, and X.q V 
as the “centered” AVE data cube set, in which the scan-temporal- mean is removed from the initial data cube 
set. 

x c y = G(n,k,x,y) -M G(k,x,y). (1) 

Similarly, the notation for the accurately calibrated AERI radiances is given as A(n,k). Since AERI has one 
pixel element, the variables x, y are eliminated from the following expressions. 

A (*o = and ( 2 ) 

n 

X A = A(n,k)-A(k). (3) 


2.2 PC Regression Calibration 

The principle component regression model that relates the GIFTS AVE data cube set and AERI data set is 

X A = X% y R. x ’ y + F, (4) 



where H x,y denotes the GIFTS- AERI regression matrix and F represents random noise and errors. By using the 
singular value decomposition (SVD), and can be factorized into 


X A =T A Pl and X.^ = (t g PI) X ’ V . (5) 

We define the number of AERI and GIFTS eigenvectors to be apply into the prediction model as variables a and 
g , respectively. The values of a and g are set to be four, which were determined experimentally such that both 
the individual features of GIFTS pixels and AERI are preserved in the prediction model. The matrices X^ a and 
for AERI and GIFTS can be related by X^ a = (TG g C ) x,y ; the matrix C x ’ y is estimated by applying the 
least squares inversion, i.e., 

c x ’ y = ((T£ s T gJ- 1 Tg s X Aa y ,V . (6) 

Finally, we can predict the GIFTS radiances G(n, k, x, y ) using the GIFTS- AERI regression matrix = 

p£VC x ’V : 

G(n, k, x, y) = A(k) + \G(n , k, x, y ) - T?(k, x, y)] (7) 
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Figure 1. FPA mean radiance images for temporal sample Set 1 (observed from Logan UT at 5 AM on 13 September, 
2006): (a) long- wave radiances averaged over 700-715 cm -1 from GIFTS original calibration, (b) shortwave radiances 
averaged over 1750-1850 cm -1 from GIFTS original calibration, (c) long- wave radiances averaged over 700-715 cm -1 from 
GIFTS PC calibration, and (d) shortwave radiances averaged over 1750-1850 cm -1 from GIFTS PC calibration. 




The FPA images for GIFTS initial calibration and PC calibration results can be seen in Figs. l(a)-(d). 
These plots indicate that the FPA electronic readout noise has been significantly reduced after applying the PC 
calibration method. Comparisons between GIFTS calibration results and AERI are made, and the differences 
in radiance between GIFTS and AERI can be seen in Figs. 2(a)-(d). After performing the PC calibration, the 
mean and standard deviation values of the differences between GIFTS and AERI radiances have been decreased. 
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Figure 2. FPA mean radiance comparisons for temporal sample Set 1 (observed from Logan UT at 5 AM on 13 September, 
2006): (a) long- wave FPA mean radiance differences between AERI and GIFTS calibration results, (b) short-wave FPA 
mean radiance differences between AERI and GIFTS calibration results, (c) long-wave FPA standard deviation of radi- 
ance differences between AERI and GIFTS calibration results, and (d) short-wave FPA standard deviation of radiance 
differences between AERI and GIFTS calibration results. 


3. TEMPERATURE AND MOISTURE RETRIEVAL 

In this section, temperature and moisture profiles are determined using the PC-calibrated upward viewing GIFTS 
radiances that are obtained through the calibration process described in the previous section. The results 
are derived and verified using radiosonde measurements collected throughout the GIFTS 9-hour sky viewing 
measurement period on September 13, 2006. The purpose of the retrieval process was to demonstrate the 
sensitivity of the calibrated GIFTS measurements to the temporal variability of atmospheric temperature and 
moisture vertical structure. For the retrieval operation, we use a new set of data cubes that combines both the 


longwave and shortwave spectral bands. Given l = 1, . . . ,L as the short-wave spectral sample number, where 
L = 1049, <ji = a i 0 + / Act, and A a = 0.5733 cm -1 ; a i covers from 1649.48 cm -1 to 2250.33 cm -1 . Define the 
new GIFTS and AERI data sets as 

= G(n,m,x,y) — G(m,x,y) and (8) 

Xi = A{n,m) — A(m), (9) 

respectively, where m = [fc, l\ combines the long-wave and short-wave spectral sample numbers. Similar to 
previous deviations, we can estimate C x,y from and X^. The mean value of G x,y averaged over the entire 
FPA, excluding bad pixels, is computed from 
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I Qx, V ' 

z ^ 

(x,y)ev 
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where V denotes the set contains the good pixels, and Z is the total count of good pixels, 
compute 

—T 


T> x ' y 
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From which, we 


(ii) 


D x,y contains the matrix of coefficients of the AERI radiance eigenvectors that are predicted from the GIFTS 
spectra. Given radiosonde measurements of temperature and mixing ratio on the same day as functions of 
pressure as T(h,ri) and Q(h, n), where n = 1,. , . , N and N = 52. The pressure level is represented by h = 
1, . . , H and H = 31. The corresponding pressure value is given by p = po + hAp; p covers from lOOOmb to 
700mb, and A p = lOmb. Define T(h) = n), T = T(h,n) — T(h) and Q(h) = jj^ n Q(h,n),Q = 

Q(h,n) — Q(h), we can establish the regression models for temperature and moisture as T x,y = (DA) x,y + Ft 
and = (DB) x,y + F q. The matrices A x,y and B x,?/ are solved using the least squares inversions: 


A x ’ y = ((D T D)- 1 D T T^’ y (12) 

W’ v = ^(D T D) _1 D T Q^ a: ' y . (13) 


The temperature and moisture as functions of pressure and time are predicted from 


T x ’ y = (DA Y’ V +T(h) (14) 

Q x ’ y = (t>bY’ V + Q(h), (15) 


respectively. Figs. 3(a)-(f) contain time cross-sections, shown as false color images, of temperature and moisture 
profiles observed by 1.5 hour interval radiosonde measurements, and the 10-minute interval AERI and GIFTS 
PC-based algorithm retrievals. The 52 10-minute time interval samples cover a time period of nine hours between 
5 AM and 2 PM local mountain daylight time at Logan Utah on September 13, 2006. In Figs. 4(a) and 4(b), the 
linearly interpolated point measurements by the radiosonde during this time period are compared to the AERI 
and GIFTS retrieved point values for temperature and relative humidity, respectively, at five different atmospheric 
levels. Comparisons are made based on these three sets of profiles; the results are shown in Figs. 4(a) and 4(b). 
The very high sensitivity of the GIFTS observations to spatial and temporal variability of temperature and 
moisture is clearly shown. The results provide confidence that if the GIFTS were to be flown on a geostationary 
satellite, the remote sensing and forecast application objectives of the GIFTS measurements would be achieved. 


4. LUNAR SURFACE TEMPERATURE AND EMISSIVITY ESTIMATION 

In this section, we discuss the results obtained from the lunar measurements collected during the moon viewing 
events (MVE). This phase of the GBM experiment is designed to capture high quality lunar images across the 
GIFTS EDU spectral channels at a scan resolution of 9.18 cm -1 . Each MVE data set is consist of multiple scans 
of the lunar images, which are radiometrically calibrated to produce the MVE radiance spectra. The atmospheric 
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Figure 3. Vertical cross-sections of temperature and moisture observed at Logan UT between 5 AM and 2 PM MDT on 13 
September 2006: (a) radiosonde temperatures, (b) AERI predicted temperatures, (c) GIFTS PC predicted temperatures, 
(d) radiosonde relative humidity measurements, (e) AERI predicted relative humidity, and (f) GIFT PC predicted relative 
humidity. 
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Figure 4. Comparisons of GIFTS and AERI 10-minute interval retrieved temperature and moisture values to linearly 
interpolated 1.5 hr interval radiosonde measurements observed at five different atmospheric levels at Logan UT between 
5 AM and 2 PM MDT on 13 September 2006: (a) radiosonde temperature measurements vs. GIFTS PC predicted and 
AERI predicted temperatures, and (b) radiosonde moisture measurements vs. GIFTS PC predicted and AERI predicted 
relative humidity. 




background effects are then removed from scans that contain lunar images. In the following step, we align the 
longwave and shortwave scans such that lunar images for simultaneous longwave and shortwave scans are spatially 
co-registering. Figs. 5(a) and (b) show the MVE radiance spectrum for a single pixel that is located within the 
lunar image and the mean brightness temperature averaged over 1100.7986-1119.1453 cm -1 , respectively. Using 
the line-by-line calculated atmospheric transmittance, we estimate the lunar surface temperature and emissivity. 
The results of estimated mean lunar surface temperature and emissivity over 1100.7986-1119.1453 cm -1 can be 
seen in Figs. 6(a) and (b), respectively. 
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Figure 5. Longwave and shortwave MVE radiance spectrum and brightness temperature: (a) MVE radiance spectrum for 
Pixel (50,70), and (b) MVE mean brightness temperature averaged over 1100.7986-1119.1453 cm -1 . 
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Figure 6. Estimated lunar surface temperature and emissivity: (a) estimated MVE lunar surface temperature, and (b) 
estimated MVE emissivity averaged over 1100.7986-1119.1453 cm -1 . 


5. CONCLUSION 

This paper describes a calibration refinement technique for the GIFTS-EDU designed to remove non-linearity 
distortions, found difficult to remove in the initial calibration of the data. This technique utilizes Principle Com- 
ponent (PC) analysis to compensate for artifacts remaining after the initial radiometric calibration process, thus 
further enhancing the absolute calibration accuracy. The algorithm is applied the data collected with GIFTS, 






together with simultaneous observations by the accurately calibrated AERI (Atmospheric Emitted Radiance In- 
terferometer), during an atmospheric measurement experiment conducted at Logan Utah on September 13, 2006. 
The algorithmic components of this method and the calibration results are discussed in the paper. Temperature 
and moisture profiles based on retrievals using the PC-calibrated GIFTS radiances are also presented. The very 
high sensitivity of the GIFTS observations to spatial and temporal variability of temperature and moisture is 
shown. Furthermore, using the GIFTS GBM calibration model, we compute the calibrated radiances from data 
collected during the moon tracking and viewing experiment events. From which, we derive the lunar surface 
temperature and emissivity associated with the moon viewing measurements. The results provide confidence 
that if the GIFTS were to be flown on a geostationary satellite, the remote sensing and forecast application 
objectives of the GIFTS measurements would be achieved. 



